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' New solutions for second-order intertwining relations in two-dimensional 

, SUSY QM are found via the repeated use of the first order supersymmet- 
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rical transformations with intermediate constant unitary rotation. Poten- 
tials obtained by this method - two-dimensional generalized Poschl- Teller 
potentials - appear to be shape-invariant. The recently proposed method 



. of 5'[/5y— separation of variables is implemented to obtain a part of their 



spectra, including the ground state. Explicit expressions for energy eigenval- 
ues and corresponding normalizable eigenfunctions are given in an analytic 
form. Intertwining relations of higher orders are discussed. 
PACS numbers: 03.65.-w, 03.65.Fd, ll.30.Pb 

1. Introduction 

The importance of each new exactly solvable model in one-dimensional (ID) Quantum 

Mechanics is well known, especially because the list of such models is quite small. The 

elegant modern approach used for the study and classification of these "elite" models was 

provided by Supersymmetrical Quantum Mechanics (SUSY QM) PP,|2I, which is in essence an 

alternative formulation of the famous Factorization Method jS| in one- dimensional Quantum 

Mechanics. Furthermore, the introduction in the framework of SUSY QM of a new notion - 

the shape invariance - gave a novel, algebraic, tool to deal with such kind of models. 

There are different ways of going beyond the scope of the standard Witten's SUSY QM 

in order to enlarge the class of involved models. The Higher Order SUSY QM (HSUSY 

QM), or, equivalently. Polynomial and iV— fold SUSY QM |S],|n! as well as constructions for 

multidimensional coordinate spaces [Zj,|H] are among the most promising ones. 
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From the very beginning, after ID SUSY QM was formulated by Witten the question 
of finding the opportunity to generahze it for higher dimensions of space attracted consider- 
able attention. A direct rf— dimensional generalization was built in ^ by means of methods 
originating from SUSY Quantum Field Theory. In this approach the Superhamiltonian (of 
block-diagonal form) includes both scalar and matrix components and can be used to analyse 
different physical problems with matrix potentials 0. 

In the particular case of two-dimensional space an alternative SUSY QM approach was 
proposed, which directly generahzes the HSUSY QM ideas, namely, the use of the SUSY 
intertwining relations with the second order supercharges. This method avoids an appearance 
of matrix potentials and provides the intertwining of two scalar Schrodinger Hamiltonians. 
A large class of such intertwined Hamiltonians was found in jTHj - [T!^ . 

In the framework of the latter approach two new methods for the study of the spectra and 
the (normalizable) eigenfunctions of the two-dimensional quantum models were proposed re- 
cently in ^3], ^n]: 5't/S'y— separation of variables and the two-dimensional shape invariance 
(see also the review- like paper ) • The combination of both of them was explored to inves- 
tigate a specific model - a generalized 2D Morse potential with three free parameters - which 
is not amenable to the conventional separation of variables. As a result, this model turned 
out to be partially solvable, i.e. only a part of the variety of its normalizable wave func- 
tions and corresponding values of energies were found analytically. Thus the transfer from 
one-dimensional to two-dimensional shape invariance was accompanied by the loss of the 
complete solvability with only the partial one remaining. It is worth mentioning here that 
each of the 2D Hamiltonians involved in the second order intertwining relation is integrable: 
the symmetry operator of the fourth order in derivatives was constructed explicitly in terms 
of supercharges jTU], [12] • 

In this paper both approaches of two-dimensional SUSY QM - the direct two-dimensional 
generalization |Z|,jH] and the second order construction [ID]- IE] - will be used to build and 
to investigate some new models, to which no standard separation of variables can be applied. 
Again Sf/ 5*1^— separation of variables turns out to be applicable to the model, providing a 
set of normalizable wave functions. This model, which is shown to be partially solvable, will 
be called a 2D-generalized Poschl- Teller potential. 

As for the method of the two-dimensional shape invariance - the situation is 
more delicate. Though the considered model possesses the property of shape invariance, the 
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corresponding solutions of the Schrodinger equation turned out to be unnormalizable. 

The paper is organized as follows. In Section 2 the known methods of 2D SUSY QM 
will be described briefly in order to simplify the comprehension of the new results. A new 
technique of searching for solutions of the two-dimensional second order intertwining relations 
will be presented in Section 3, and, in particular, two-dimensional generalizations of Poschl- 
Teller potentials will be constructed. In Section 4 St/S*!^— separation of variables will be used 
to find a part of the spectrum of this model and analytical expressions for its wave functions, 
including the ground state. The peculiarities of shape invariance are also investigated. In 
Section 5 an additional structure with two different superpartners for the same Hamiltonian is 
presented, new intertwining relations of fourth and sixth orders in derivatives are constructed 
(the last ones are shape-invariant). 

2. Basics of 2D SUSY QM 



2.1. 2D representation of SUSY algebra 

The SUSY algebra of quantum mechanics is given by the following (anti) commutation rela- 
tions [l]: 

{g+,g-} = ^; {Q+,Q+} = {Q-,Q-} = 0; [Q^,H]=0. (1) 

In the case of two dimensions it can be realized [Zj, [Hj by the following 4x4 matrix operators: 
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(2) 

where two scalar Hamiltonians H'^^\H^'^^ and one 2x2 matrix Hamiltonian of 
Schrodinger type can be expressed in a quasifactorized form (compare to the factorized 
form in one-dimensional case P, 



PTPi 



-df + v'-'\x) = -df+{da{^. 



(3) 



(1) 



ik 



Qi (it +Pi Pt = -^ikdf + Sik(^{dix{x)Y - <9fx(^)) + 2didkx{x), 



with components of supercharges of first order in derivatives: 

= Tdi + dix{x); p^ = eikq^, (4) 

where di = d/dxi and summation over repeated indices is imphed. Anticommutators in 
produce the following intertwining relations for the component Hamiltonians H^^\ H^"^^ 
of the Superhamiltonian 

H^<ll = ^H^^ H^<h=<l7H^h H^p,=prH^; H^j>t=piH^. (5) 

They connect the spectrum of the matrix hamiltonian with spectra of two scalar ones. In 
general, and if^^) 

are not isospectral since q'^Pj^ = due to Q. 
2.2. Second-order supercharges in 2D SUSY QM 

Two-dimensional SUSY QM models without any matrix Hamiltonians were constructed jlUj . 
[TT] . [in] by means of second order supercharges 

= {Q-)^ = g^kimdk + C,{x)d, + B{x) (6) 

where Qik-, Ci, B are arbitrary real functions. Some particular solutions for two scalar Hamil- 
tonians which satisfy the intertwining relations 

/f(i)(x)g+ = g+ijW(f); H^''\x)Q- = Q-H^^\x), (7) 

were found. They both possess the symmetry operators R^^'"^^ of fourth order in derivatives 

cni, m 

which are not, in general, polynomials of H^'-\ 

In terms of unknown functions gik,Ci, B,V^^\V^^^ Eg.© has the form ^0] of seven 
nonlinear partial differential equations, and its general solution is not known. To obtain 
particular solutions different ansatze for "metrics" Qik were used. 

Only the choice of Laplacian (elliptic) metrics gik{x) = diag{l, 1) leads to Hamiltonians 
amenable to i?— separation ^\ of variables. All other choices of metrics give nontrivial 
results. The case of Lorentz (hyperbolic) metrics gik = diag{l, —1) was investigated in 



4 



papers (TU] - [IHl- In particular, the intertwining relations (|7j) were reduced to the pair of 
differential equations: 

91^2^ = 0; d4C-F) = -d+{C+F). (8) 

where'^ x± = (xi ±^2) / and Ci^2 were proven to satisfy C± = Ci^Ca = C±{V2x±). Then 
potentials and the supercharges are expressed in terms of functions C±{^/2x±) 

and F{x) which obviously can be written as F = Fi{2xi) + ^2(2x2) according to (jH)): 

+^(f2(2x2)-Fi(2xi)); (9) 
Q+ = (dl-di) + Cidi + C2d2 + B; (10) 
5 = i(c+(v^a;+)C_(v^x_)+Fi(2a;i) + F2(2x2)), (11) 

where the prime denotes the derivative of function with respect to its argument. A list of 
particular solutions of (jH)) was obtained in pi] - [12] . In the next Section we will obtain new 
solutions for the case of hyperbolic metrics. 



3. New solutions for the Lorentz (hyperbolic) metrics 



3.1. Intertwining of second order with reducible supercharges 

Let us consider two Superhamiltonians H and H of 2D SUSY QM: 
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( i7(°)(f) \ 

ffi;)(x) 

\ H^'^\x)] I 
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H^''\x)- j 



(12) 



with superpotentials xi.^) ^-nd x{.^)i correspondingly. 



In addition, let and HH.' be linked by an unitary 2x2 matrix transformation U: 



r(l) 



TT S-(l) _ TtWtt . 



U = aoao + ia a ; + a =1; ao, ctj ^ 



(13) 
(14) 



■^We use here the definition of x± sHghtly different from the analogous one in JT] - |13| . 
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where cxj are the Pauh matrices and (Tq is the unit matrix. 

Then (due to (0)) the scalar Hamiltonians H^'^^ and H^^^ can be included in the chain: 

i7(o) 4U ijW H^l^ A #(0), (15) 
leading to the intertwining relations between a pair of scalar Hamiltonians: 

ifWg- = g-^w, Q+i/W = ijWg+ (16) 

with second order operators 

Q- = (g+)t = g+f/,,g-. (17) 

This intertwining operator is constructed from two first order ones with intermediate matrix 
transformation Uik- Precisely this matrix provides that such supercharges are nontrivial 
and, contrary to Subsection 2.1., can be naturally described as reducible (compare with 
the case of one-dimensional reducibility introduced in ^5j). 

In contrast to the approach of [13] (see Subsection 2.2.), the first Hamiltonian in the 
chain (fT3|) is quasifactorized according to Eq.©. Therefore the solution of the corresponding 
Schrodinger equation with zero energy can be written as "^^^^ ~ exp (— x)- Due to expression 
fll7|) . exp (— x) is a zero mode of supercharge as well. In general, until the specific form 
of x(^) is chosen, the normalizability of this solution is not guaranteed. But in the concrete 
model (j3T|l analyzed below in Subsection 3.3. the zero energy solution is normalizable due 
to asymptotic properties of x{^) corresponding ranges ()4(i|) of parameters. 

Target Hamiltonians H^^^ and H^^^ are expressed in terms of two unknown functions x 
and X (see ©)• To determine these functions one should substitute ©, (HH) and (fT7|) into 
fll(i|l . After some manipulations one obtains the system of equations for x± = (x i x)/2: 

ttaDX- + 2aidid2X~ = 0; aiOx+ - 2a-idid2X+ = 0; (18) 
a2Clx+ + 2aodid2X- = 0; aoDx- + 2a2did2X+ = 0; (19) 
{dkX-){dkX+) = 0, (20) 

where □ = dl — dl- Eq. fEIH) (which is equivalent to {dkxY = {^kxY) can be used to simplify 
expressions (0) for the Hamiltonians if*-"-* and H^^^: 

^(0)^^(0) ^ -d^ + {{da^f - dh+) + {{da-fTdh-). (21) 

Linear partial differential equations (fTH)) -(fT ^ can be easily solved, but the solution of the 
nonlinear Eq. ()20j) is a nontrivial problem. 
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3.2. The particular solutions of the intertwining relations 

It can be shown that in the case when all coefRcients a, and ao in (|T4|l do not vanish, 
potentials V^^^ and V^'^^-' are 4-th order polynomials on xi^2 with some additional constraints 
for their coefficients. In the present paper we will consider potentials beyond this rather 
narrow class, restricting ourselves to the particular case ao = ai = ^2 = 0; 0:3 7^ 0, 
i.e. U = (73. Then the metrics of supercharges is Lorentz, i.e. they belong to the class 
discussed in Subsection 2.2. 

For this case equations (fTSj) - (fT^ read: 

Dx- = 0; d,d2X+ = 0. 

Their solution is 

X+ = /ii(a;i) +yU2(a;2), 
with /ii^2, At± being arbitrary functions. The last equation (jSDI) takes the form 

n[{xi) +/U2(a;2) [f^'+ix^) - = 0. 

By substitutions = /i', it becomes purely functional (without derivatives) equation: 

+ 0_(x_)] = -Mx2) [4>+M - • (22) 

The general solution of is given in the Appendix'^. Some particular cases will be 
discussed in Subsection 3.3. 

The Hamiltonians (PT|) and intertwining operators (fT7j) can be expressed in terms of 0: 

+ (02_(x_)T0'-(a^-)), (23) 

g± = 52 - di ± V^(0+(X+) + 0-(x_))9i T V^(0+(X+) - 0-(x_))92 - 

-(02(xi) -0;(a;i)) + (02(x2) -0'2(a;2)) +20+(a;+)0_(a;_). 

By rearrangement of terms Eq.(j22I) can be rewritten as: 

0+(x+)[0i(a;i) + 2(x2)] = -0_(a;_)[0i(a;i) - 02(0:2)], (24) 

"^It was derived by D.N. Nishnianidze (private communication). 
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i.e. in the form similar to the initial Eq.()22j). This means that ()22|) possesses the 
symmetry property (which will be called 5*1 symmetry in the subsequent text): if 
|0i(a;i), 02(3^2), 0+(a;+), 0-(a;_) I is a solution, then |0+(a;i), 0_(s2), 0i(a:;+), 02(a;_)| is 
also a solution. Let us mention one more discrete symmetry of ()2()|1 . S2 symme- 
try: {0i(a;i), 02(3:2), 0+(a;+),0-(x_) I — > { 0i(a;i), -02(x2), 0+^(a;+), 0I^(x_) I. The ^i- 



symmetry produces another supersymmetrical model: 

V(°), = (02(^+) T 0U^+)) + (02(^-) T <P'2ix-)) + [<Plix^) - <P'4x^)) + 

+ (4>1{X2)-^UX2)), (25) 

Q± = df - dl ± y2((/)i(x+) + Mx~))di T v^(0i(a;+) - Mx^))d2 - 

-(</)^(a;i) -0;(xi)) + (</)2 (0:2) -C(x2)) +20i(x+)</)2(x„). 

Below both forms (j23|l and p5|l will be explored. 

To compare the new notations of this Section with those of [TT] - [12] (see Subsection 

2.2. ) one can use the following relations: 

C±(%/2x±) = 2v^0±(x±); 

/ \ (2d) 

i^i,2(2a;i,2) = T4(0?,2(a^i,2)-0'i,2(a^i,2)). 

3.3. Nonperiodical solutions for potentials V^^\ V^^l 



From Eq. ()A8|) one can conclude that for an arbitrary choice of the parameters a,b,c the 
functions 0i^2 are expressed in terms of elliptic (Jacobi or Weierstrass) functions [18^ (volume 
3). In this paper we restrict ourselves by considering the limiting cases, for which the 
potentials are not periodical (the models with periodicity properties in xi 2 will be studied 
elsewhere). 

The integral in the r.h.s. in ()A8|) is an elementary function only if either some of coeffi- 
cients a, b, c are zero or the quadratic polynomial is a full square. There are two families of 
solutions of ()22j). with members interconnected by symmetries Si and 5*2 (and their combi- 
nations) : 

(i) (pi{x) = (t>2{x) = A/x] 0+ = 0_ = B/x {A,B = const). 

The multiparticle potentials of this type were found in to be quasi-exactly solvable j2Uj . 
All other members of this family allow the separation of variables. 

(ii) 01 (x) = 02 (x) = M + 5_e""") ; 



-L 



For the particular case 5_ = one has: 



= Lcoth ax/V2 . (27) 



+ A(2A T a) [sinh (f (xi - Xa))] 



(2^ 



with two new constants A, B instead of M, L, 5+. This potentials (up to translations in Xi^2) 
were analyzed in ^3] and were found to be shape-invariant. 

Another particular case 5+ = —5- for ()27j) after using symmetries and redefinition of 
parameters gives: 

sinh \/2ax 

0+ =(/)_= Btanhax, (29) 

The corresponding potentials and intertwining operators for this model due to (j23|l are: 

B(B±a) \ / o B(B±a) 



y(o)^y(o) 



B' 



+A 



cosh 



■V2 



{xi + Xa)) 



+ B' 



cosh (-^(xi - X2)) 



+ 



A — \/2a cosh {\/2axi) _|_ ^ + V^ctcosh (v^axa) 



sinh^ (v/2q;xi 
- 9^ ± V2fi 

tV25 

y4 — cosh (-\/2«a;i 



-A 



sinh^ {\/2ax2) 
tanh ^— ^(xi + ^2)^ + tanh ^— ^(xi — X2) 

v4 + v^a cosh {^/2ax2) 



(30) 



5iT 



sinh^ (V2 



sinh^ (-\/2q;X2) 



+ 



/ a 



+2B tanh ^— ^(xi + X2)J tanh ^(xi — X2 
Solution obtained by the discrete symmetry Si, is: 



52 



+v4 



fi(fi + a) 
cosh (axi] 



B' 



B{B + a) 
cosh^ (0x2) 



+ 



A =1= v^a cosh (a(xi + X2)) _|_ ^ v^*^ cosh. {a{xi — X2)) 



Q± = dl-dl±V2A 



sinh (q;(xi + X2)) 
1 



+ 



sinh (a(xi — X2)) 
1 



(31) 



sinh (a(xi + X2)) sinh (a(xi — X2)) 



diT 
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sinh (q;(xi + X2)) sinh (q;(xi — X2)) 
B{B + a) 



cosh (axi) _ 
2^2 



+ 



B' 



B{B + a) 
cosh (ax2] 



+ 



sinh {a{xi + X2)) sinh {a{xi — X2)) 

Both potentials (jHUj) and (jHH) can be treated as superpositions of two one- dimensional 
Poschl- Teller terms plus a singular term (so we will refer to them as 2D-generalized 
Poschl- Teller potentials). Each of them possesses a term which prevents application 
of the conventional method of separation of variables to determine their eigenfunctions and 
eigenvalues. Meanwhile, a part of the spectrum and corresponding eigenfunctions will be 
found by the method of SUSY— separation of variables (see fl], ^Hj) in the next Section*^. 



4. Partial solvability of 2D-generalized Poschl- Teller po- 
tentials 

4.1. SU SY —separation of variables 

As far as Hamiltonians with potentials (p?T|) are intertwined by operators with Lorentz 
metrics, we shall briefly remind the reader of the general method for searching for eigenvalues 
and eigenfunctions proposed in ^3]. From intertwining relations Q'^Ti.^^^ = Ti.^^'^Q^ (where 
7-^(0) = — _|_ y(o) g^j^^ ^{0) ^ _q2 _|_ y(o)-) obtains that the subspace of zero modes^ of 
the supercharge 

Q+Q{x) = 0, (32) 

is closed under the action of H^^^: 

H^°^Q{x) = CQ{x) (33) 

with some constant matrix C. 

To determine the eigenvalues Ek and eigenfunctions \E'fc(5) of H^^^ one needs (see more 

details in jHj) a matrix B, which satisfies the matrix equation BC = KB with an unknown 

^Other members of the same family {S2- and 5*2 S*! -symmetric to ((23) can be treated analogously. 
^Here we suppose that {N + 1) normalizable zero modes r2„(a;) are known, and il(x) is a column vector 
with components f2„(ai'), n = 0, 1, ... TV. 
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yet diagonal matrix A = diag{Xo, Xi, ■ ■ ■ , Xn)- Then the matrix B transforms zero-modes 
QnS into wave functions \l/„'s. 

Operator belongs to type For this type of supercharges the problem (IH2|1 permits 
the conventional separation of variables in by means of "gauge" transformation, which 
separates variables in the supercharge: 



M^)Q+^<x) ^ q2 _ g2 ^ l^F,{2Xi) + i"2(2X2)), 

1 



(34) 



h{x) = e-"("^)7^(°)(^)e"^"^ = -d! - dl + - C2{x)d2 - -tF,{2xi) + -^2(2x2). (35) 



j C+{V2x+)dx+ + J C-{V2x-)dx. 
uJx) = e-'^^^^aJx). 



Then the zero modes u!n{x) of q+ can be written as products co'„(x) = rin{xi)pn{x2), where 
Pn and T]n are eigenfunctions of the one-dimensional Schrodinger equations with "potentials" 
(=F|-Fi,2(2a;i,2)), correspondingly (see (jBll)), and common eigenvalues (constants of separa- 
tion) e„. 

It is obvious that huj = Cuj with the same matrix C as in (jHHj) . The simplest way to find 
C is to calculate the r.h.s. of (jHHjl . which can be rewritten as: 



huj„ 



[2e„ + Ci(f)9i-C2(x)92]cu„. 



(36) 



As a result, after construction of the matrix B one will obtain part of the spectrum Ej, and 
corresponding wave functions \E'fc(a;). 



4.1.1. Calculation of C 

This general method, proposed in [T3] and used there successfully to investigate the 2D Morse 
potential, can be applied to the pair V^^-*, V'-''-* as well. In this case both one-dimensional 
equations for multipliers rin{xi) and p„(x2) have the same "potentials" - one dimensional 
Poschl- Teller potentials - being exactly-solvable: 



-dl + B'^ - B{B + a) cosh ^ (axi))?7„(xi) = e„77„(xi) 



(37) 



and a similar equation for Pn{x2)- By the change of variable ^ = tanhaxi Eq. ()37|1 can be 
reduced to the generalized Legendre equation |21j : 



d 



dr] 



+ 
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where s = B/a. To have the finite solution at ^ = — 1 the condition a/ (-B^ — e)/a'^ — s = —n; 
n G N must be satisfied. It gives a discrete set of values for the separation constant e: 

en = a^n(2s — n) 

for n < s. Corresponding functions (up to normalization factors) are rin = Pg~^{^), where 
Py{x) are the (generalized) Legendre functions. Thus, one achieves the expression for ujn'- 

^n = i^r"(6)pr"(6), 

with = tanhaa;^, i = 1,2. 

The next step in calculating the eigenfunctions is evaluating the r.h.s. of (jHUjl : 

(I _ p2\l/2Q _ f2\l/2 

hunix) = 2a^n{2s-n)uJn-2V2Aa{2s-n){n+l)^ ^f-^^-^f^^U{n, s; ^i, ^2) , (38) 

where me made the shorthand notation 

n(n,s;ei,6) = 6(e? - iy^'pr'''\^i)pr'i^2)-ue2 - i)^/'pr"(ei)pr"-^(6)- 

Our goal is to represent the r.h.s. of (|38p as a linear combination of uj^s. For this purpose 
we use the reccurrent formula for Legendre functions (see ^H], v.l, p. 161, eq. (1)) 

Applying it twice to n(n, s; ^1; .^2) one obtains the following reccurrent formula: 

1 



n(n,s;^i,^2) 



{n + l)n{2s - n){2s - n + 1) 



2is~n + l)ie,-e) 

(1_ ^2)1/2(1 _^2)i> »-i+n(^-2, s- 6, 6) 

(39) 



To stop this procedure at n = 0, one has to consider s G N. In this case the Legendre 
functions turns into the associate Legendre polynomials, for which P^{z) = for m > n. 
So, applying several times. 



n(-^ ^' - (i-ef)v^(i-el)v^ S 



fc=0 



with constants a„fc. The matrix elements Cnk of matrix C are: 

Cnfc = 2a^n(2s - n)5nk - 2\f2Aa{2s - n) (n + l)a„fc; (40) 

0, k > n; 

ank={ 0, A; = n-2m-2; m = 0, 1,2, ... (41) 

(fc-l)!(2s-n-l)! 



r,/ _ J \ (K-L)'.(ZS-n-i 

^{^ (n+l)!(2s-fc)! ' 



k = n — 2m — 1; m = 0, 1, 2, 
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4.1.2. Calculation of eigenfunctions 



Matrix C for the model ()3H1 appeared to be triangular, and hence its eigenvalues coincide 
with the diagonal elements: 

Ek = Ckk = 2a''k{2s-k). (42) 

This formula gives Eq = 0, demonstrating that the zero energy solution of TC^^^ is a zero 
mode of Q'^ as well: \t'o ~ exp {—x{x)). 

In order to avoid zeros on the diagonal of C one can shift Hamiltonians by a constant 
7. This transformation does not destroy the intertwining relations (|T6|l and changes C as 
follows: Cik Cik + 75jfc- This new C can be diagonalized by the method, presented in fT3]. 
Namely, the formal solution for the matrix elements of B reads: 



'N-p-1 

E 

. 1=1 



(43) 

N~m,p 



where (A^ + 1) matrices t^™'^ are defined by 



(m) Cn^k 



Tji k — 

CN-m,N-m ~ Ck,k 

and label (m) has values m = 0,1,. ..,A^. In Eq. fj43|) the repeated index — m is not 
summed over, and (to avoid misunderstanding) r*^™) means the /th power of matrix r*^™'-'. 
Thus one obtains the recipe for the construction of eigenfunctions for 7^(0) in (jSD): 

N 

'^N-ni.x) = bnAix). (44) 

Formula gives us the opportunity to express an element bm,p by means of the r^™) 
matrices and an arbitrary element bm,N-m on the crossed diagonal. This last element can be 
fixed by the normalization condition for "^N^m- The reason for the "inverted" numeration 
of \I' in ()44|) is to make "^k dependent only on Qf, I = 0,1,. ..,k. In particular, \E'o ~ fig. 
So, applying the method of St/ S'F— separation of variables to H^^\ H^^^ one obtains a set 
of eigenvalues and eigenfunctions for 'H^'^\ 

Keeping in mind that s = B /a > 0, we can restrict ourselves with B > 0, a > 0. The 
conditions of normalizability of Qn (and therefore, of ) for all oOn can be derived from the 
explicit expressions: 

f (cosh {a{xi + X2)) - l)(sinh(tt(xi -xa))) ^ ^/(v^a) ps^uf. x 

(]„=a;„exp«:= ^ sinh (a(xi + x,))(cosh (a(xi - x,)) - 1) ) ^^'^^^ 

(45) 
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The constraint is: < „ < The condition A > keeps the strength of attractive 

singularities of both superpartners V'-"-', V^^-' at x± — not exceeding the standard bound 
— l/(4x^). The resulting range of parameters for both Qn and fin is: 

B a 
a>0; B>0; — G N; < A < —. (46) 

a V2 

At first it seems that the energy eigenvalues ()42j] . which were built above from the analysis 
of the zero-modes of Q^, should be absent in the spectrum of its superpartner 'H^^\ since 
the corresponding eigenfunctions are annihilated by Q^. However, the whole procedure of 
Sf/SF— separation of variables can also be implemented for the spectral problem for Ti.^^^ 
by suitably replacing ()H2j) with Q^Q = 0. Since Q~ and Q'^ differ only by sign in front of 
the first derivatives (see ()25p). one should use the "gauge" transformation with exp (— k(x)). 
In this case one will obtain the same equations ()37p. as for the problem ((221) • Then the 
zero-modes of Q~ can be written as: 

fln{x) = exp { — K{x))LJn{x) = CXp ( — 2/«(x) (^) , 

and the corresponding matrix C is again triangular. To be more precise, it is the same as 
fl4U|) - (PT|) up to the sign of the last term in pn|) . Therefore, its eigenvalues, i.e. values of 
energy for 7i^°\ coincide with (021) • One can check that the eigenfunctions are normalizable 
in the same range of parameters ()4(ij) . Thus the obtained part of the spectra of superpartners 
H^'^^ and H.^^^ totally coincide. In a certain sense this result is similar to one of the variants 
of the second order intertwining in ID HSUSY QM [S]: the equal number of bosonic and 
fermionic zero modes does not signal the spontaneous breaking of the supersymmetry. 



4.2. The method of shape- invar iance 

Shape-invariance j4j , [2j , pjy is an additional property of intertwined superpartner Hamilto- 
nians which gives the opportunity to determine their spectra algebraically. Namely, if both 
Hamiltonians depend on some extra parameter (or set of parameters) a, this property reads: 

H{ao) = H{ai)+n{ao), (47) 

where ai = /(oq) is another value of parameter, and 7l{ao) does not depend on x, i.e. H 
has the same (up to an additive constant) shape as H, but with another set of parameters. 
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Let us assume that we know some eigenfunction "^^^^ of H and the corresponding eigen- 
value in some range of the parameter a. Starting from 

H{ai)¥°\ai) = EW(ai)^(°)(ai) (48) 

and employing (jTfjl . one obtains: 

H{ao)¥'\a,) = (E(°)(ai) +7^(ao))vI/(°)(aO. (49) 

Using intertwining relations (fTBj) in (j^ . 

H{ao) [Q-{ao)¥'\ar)] = (^^(ai) + 7^(ao)) [Q- {ao)¥'\ar)] . (50) 

This means that H{ao) has the eigenvalue E^^\aQ) = E^'^\ai) +7^(ao) with the wave function 
^^^■'('^o) = Q" {0'o)^^^Kcbi) (its normalizability is not guaranteed). Starting from PH|) with 
parameter a2 = /(/(ao)) and repeating the described procedure twice one can find: 

H{ao) [Q-{ao)Q~{ai)¥''\a2)] = (E(°)(a2)+7^(al) + 7^(ao)) [Q' {ao)Q- {a,)¥'\a2)] , (51) 

which gives one more point (\E'''^^(ao), -E'^^-'(ao)) in the spectrum of H^ao). The general for- 
mulas are 

^^"^ao) = Q-{ao)Q-{a^)...Q-{a^^^)¥'\an), (52) 

n-l 

E(")(ao) = i?(°)(0 + E^(«^)- (53) 

k=0 

So, we have constructed a " shape- invariance chain" of eigenf unctions starting from 
one given. The natural idea is to combine this method with 5* f/ 5*1^— separation of vari- 
ables (if the Hamiltonians possess shape-invariance): having (N+1) eigenfunction from 
S'f/S'F— separation, we use each of them to start the described above shape-invariance chain. 
This procedure was implemented for the generalized 2D Morse potential ()28|1 in ^3]. 

For the 2D-generalized Poschl- Teller potential the situation becomes more complicate. 
Indeed, the Hamiltonians H^^\ H^'^^ with potentials (jHUj) are shape- invariant: 

(f ; B, a) = {x;B-a,a) + 2{B^ -{B-af)), (54) 

where, in the notations introduced above, Oq = B; ai = /(oq) = B — a. But, contrary to 
the model ^3], the method of S'f/S'y— separation of variables does not work here, since the 
zero modes Un are unnormalizable for all values of the parameter A. 
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This obstacle can be overcome by exploring the relation between systems fjSUj) and (jSD): 

if(o)(xi,x2) = n^'^\x+,x.), (55) 

where in the r.h.s. arguments xi,X2 are substituted by Because a part of the 

spectrum (and the eigenfunctions) of was found by the method of Sf/^y— separation 
of variables (Subsection 4.1.), one can use the relation to obtain the corresponding part 
of the spectrum (and the eigenfunctions) for H^^\ Then one can use these data to start 
shape invariance chains for the system H^^^ , H^'^^ according to (jK^ with operators Q~ . For 
example, starting from the first zero mode one obtains: 

¥''\xi,X2;B) = Q'{xi,X2;B)...Q~{xi,X2;B - {n - l)a)Qo{x+,X-; B - na). (56) 

The general formula for the spectrum can be obtained from (j53|) . where the E^^^^s for 
each chain are taken from ()42|1 : 

Emn = 2a^[m(2s - 2n - m) + n{2s - n)] = 2a^{m + n)[2s - (m + n)], (57) 

where < m < s corresponds to the number of the chain (number of eigenfunction con- 
structed by Sf/SF— separation), and < n < s in order to keep positive all of 7l{ak), 
k = 0, (A^ — 1), since the ground state energies E(o)(afc) = 0. Comparing (jSH) with ^ 
one will find that these points of the spectrum coincide exactly {k = m + n). But at closer 
examination all "wave functions" of the form with n > 1 are unnormalizable due to the 
singular behaviour of the supercharges (j3n|l at xi,2 — ^ 0. Thus the seeming {k + 1)— fold de- 
generacy of A;— th energy level in (jHTj) is spurious since only one of the solutions of Schrodinger 
equation (namely, the linear combination of zero modes f2„) is normalizable. Therefore, in 
contrast to the method of SU SF— separation of variables, the method of 2D shape invariance 
is powerless to give normalizable shape invariance chains of wave functions for the 2D 
Poschl- Teller potential. 

5. Two-dimensional intertwining relations of more than 
second order 

In this Section we will imply equivalence of H^'^^ and 7^'-°^ up to a change of variables for the 
new construction. Due to ()55p . the intertwining relation Ti.^'^^Q^ = Q^'H^'^^ can be rewritten 
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as: 



where Q^(xi,X2) = and 7ll'^°)(xi, X2) = 7i^°-'(x+, a;„). Comparing it with inter- 

twining relations for the pair one can conclude, that H^^^ has two different 

superpartners: 



7^(0) 



(58) 



intertwined by different supercharges. Therefore, the Hamiltonians H^'^^ and i/^'^^ can be 
considered as superpartners intertwined by the fourth order operators Q^Q^. This pair 
does not obey the shape-invariance property. 

Because, the Hamiltonian H^^'' is shape- invariant ()47|1 . one can develop the construction 



= (ao) 



Hi^)^ao) H^'\ao) = H^'\a^) + n{a,) 1^'^' H^'\a{) + 7^(ao), (59) 

where ao = -B, ai = B — a (see Subsection 4.2.). The outermost operators in (jKHjl are 
intertwined by sixth order supercharges according to: 



7i(°n«o) Q+(ao)g-(ao)Q-(ai) = {a^)Q- {a^)Q- {a^) n^'\a,)+n{ao 



(60) 



and contrary to the previous, fourth order, case, obey the shape-invariance property. Thus 
one can continue the construction of the spectrum of H^^^ 



Appendix. The general solution of the functional equa- 
tion 



Applying operator {df — 9|) to both sides of (j^ . one has: 

1 (0-(^-)-0+(^+)) = 



0i(xi) 02(3^2) / V02(a;2) M^i) 

where the notation 4>i{xi) = l/0i(xi) was used. The general solution of ()A1|) is: 



-{x-) -0+(x+)), 
(Al) 



0i(xi)(0;(xi)0'2(x2))-^/^A MEilrf:,^ _ J Mgrf:,, 
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(A2) 



where A is an arbitrary function. The corresponding expression for + 0+(x+)) can 

be obtained using initial Eq.| 



Expressions for (f)± in terms of A should depend on proper argument, i.e. 9-t0=p = 0, leading 
to additional constraints for the function A: 



2 \(f)2{x2)(j)'2{x2) \4>'i{xi) 02(a;2) 

201 (Xi) 202(X2) / \MX2) MX1)J 

Its trivial solution A = gives 0+ = 0_ = 0, 0i_2 - arbitrary, and the potentials and 
(j25|l are amenable to separation of variables. 
Otherwise one can exclude A from ()A4|) - ()A5|) : 

0l(Xi)0^(x2) 02(a;2)0?(Xi) 2/ ^ \ , A / \ A ( \ / A «^ 

—J—. —J—. = 202 {X2) - 201 [xi) + 0i(a;i)0i(a;i) - 02(x2)02(a;2)- (A6) 

01 (xi) 02(a;2) 

Though there is no separation of variables in ()A6|) . it will appear after applying the operator 
did2, so that: 

m<t>^y _ m<^-y = 2a - const (A7) 

Integrating, multiplying by 0'^ 2? integrating again and taking into account ()A6|1 . one obtains 
the general solution of (|^^ in the form: 



b'l = a0i 2 + 60?,2 + c; x = ± !,9 =i ^' c = const. {A. 
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